#### 10/26/2018
#### Replication Code:
#### Goldman, Seth and Daniel J. Hopkins
#### "Past Place, Present Prejudice"
#### 
#### R code enabling analyses of
#### 2008 National Annenberg Election Studies
#### / 2012 Institute for the Study of Citizens and Politics 
#### survey

##### Load libraries
#####
library(texreg)
library(lme4)

##### load data set
load("panel1-8nhw06082017-retro-forcontextanal.Rdata")

#####
dta3$CODEMPCT10 <- dta3$DEMVOTEPCT

##### Recode stereotype data
dta3$BLKHRDWRK3 <- NA
dta3$BLKHRDWRK3[! dta3$PPETHM_3=="Black, Non-Hispanic" & ! dta3$PPETHM_3 %in% c(NA)] <-as.numeric(dta3$DETH2_A_3[! dta3$PPETHM_3=="Black, Non-Hispanic" & ! dta3$PPETHM_3 %in% c(NA)])
dta3$BLKHRDWRK3[dta3$DETH2_A_3 < 0] <- NA

dta3$BLKINTEL3 <- NA
dta3$BLKINTEL3[! dta3$PPETHM_3=="Black, Non-Hispanic" & ! dta3$PPETHM_3 %in% c(NA)] <-as.numeric(dta3$DETH2_B_3[! dta3$PPETHM_3=="Black, Non-Hispanic" & ! dta3$PPETHM_3 %in% c(NA)])
dta3$BLKINTEL3[dta3$DETH2_B_3 < 0] <- NA

dta3$BLKTRST3 <- NA
dta3$BLKTRST3[! dta3$PPETHM_3=="Black, Non-Hispanic" & ! dta3$PPETHM_3 %in% c(NA)] <-as.numeric(dta3$DETH2_C_3[! dta3$PPETHM_3=="Black, Non-Hispanic" & ! dta3$PPETHM_3 %in% c(NA)])
dta3$BLKTRST3[dta3$DETH_C_3 < 0] <- NA

dta3$BLKSCALE3 <- dta3$BLKHRDWRK3+dta3$BLKINTEL3+dta3$BLKTRST3

dta3$BLKHRDWRK4 <- NA
dta3$BLKHRDWRK4[! dta3$PPETHM_4=="Black, Non-Hispanic" & ! dta3$PPETHM_4 %in% c(NA)] <-as.numeric(dta3$DETH2_A_4[! dta3$PPETHM_4=="Black, Non-Hispanic" & ! dta3$PPETHM_4 %in% c(NA)])
dta3$BLKHRDWRK4[dta3$DETH2_A_4 < 0] <- NA

dta3$BLKINTEL4 <- NA
dta3$BLKINTEL4[! dta3$PPETHM_4=="Black, Non-Hispanic" & ! dta3$PPETHM_4 %in% c(NA)] <-as.numeric(dta3$DETH2_B_4[! dta3$PPETHM_4=="Black, Non-Hispanic" & ! dta3$PPETHM_4 %in% c(NA)])
dta3$BLKINTEL4[dta3$DETH2_B_4 < 0] <- NA

dta3$BLKTRST4 <- NA
dta3$BLKTRST4[! dta3$PPETHM_4=="Black, Non-Hispanic" & ! dta3$PPETHM_4 %in% c(NA)] <-as.numeric(dta3$DETH2_C_4[! dta3$PPETHM_4=="Black, Non-Hispanic" & ! dta3$PPETHM_4 %in% c(NA)])
dta3$BLKTRST4[dta3$DETH_C_4 < 0] <- NA

dta3$BLKSCALE4 <- dta3$BLKHRDWRK4+dta3$BLKINTEL4+dta3$BLKTRST4

dta3$BLACKGOVJOBS_1S <- NA
dta3$BLACKGOVJOBS_1S[dta3$blackgovjobs_1=="Strongly Agree"] <- 4
dta3$BLACKGOVJOBS_1S[dta3$blackgovjobs_1=="Somewhat Agree"] <- 3
dta3$BLACKGOVJOBS_1S[dta3$blackgovjobs_1=="Somewhat Disagree"] <- 2
dta3$BLACKGOVJOBS_1S[dta3$blackgovjobs_1=="Strongly Disagree"] <- 1

dta3$BLACKGOVJOBS_7S <- NA
dta3$BLACKGOVJOBS_7S[dta3$blackgovjobs_7=="strongly agree"] <- 4
dta3$BLACKGOVJOBS_7S[dta3$blackgovjobs_7=="somewhat agree"] <- 3
dta3$BLACKGOVJOBS_7S[dta3$blackgovjobs_7=="somewhat disagree"] <- 2
dta3$BLACKGOVJOBS_7S[dta3$blackgovjobs_7=="strongly disagree"] <- 1

dta3$BLACKGOVJOBS_8S <- NA
dta3$BLACKGOVJOBS_8S[dta3$blackgovjobs_8=="strongly agree"] <- 4
dta3$BLACKGOVJOBS_8S[dta3$blackgovjobs_8=="somewhat agree"] <- 3
dta3$BLACKGOVJOBS_8S[dta3$blackgovjobs_8=="somewhat disagree"] <- 2
dta3$BLACKGOVJOBS_8S[dta3$blackgovjobs_8=="strongly disagree"] <- 1


dta3$BLACKGOVSPND_1S <- NA
dta3$BLACKGOVSPND_1S[dta3$blackgovspen_1=="Strongly Agree"] <- 4
dta3$BLACKGOVSPND_1S[dta3$blackgovspen_1=="Somewhat Agree"] <- 3
dta3$BLACKGOVSPND_1S[dta3$blackgovspen_1=="Somewhat Disagree"] <- 2
dta3$BLACKGOVSPND_1S[dta3$blackgovspen_1=="Strongly Disagree"] <- 1

dta3$BLACKGOVSPND_7S <- NA
dta3$BLACKGOVSPND_7S[dta3$blackgovspen_7=="strongly agree"] <- 4
dta3$BLACKGOVSPND_7S[dta3$blackgovspen_7=="somewhat agree"] <- 3
dta3$BLACKGOVSPND_7S[dta3$blackgovspen_7=="somewhat disagree"] <- 2
dta3$BLACKGOVSPND_7S[dta3$blackgovspen_7=="strongly disagree"] <- 1

dta3$BLACKGOVSPND_8S <- NA
dta3$BLACKGOVSPND_8S[dta3$blackgovspen_8=="strongly agree"] <- 4
dta3$BLACKGOVSPND_8S[dta3$blackgovspen_8=="somewhat agree"] <- 3
dta3$BLACKGOVSPND_8S[dta3$blackgovspen_8=="somewhat disagree"] <- 2
dta3$BLACKGOVSPND_8S[dta3$blackgovspen_8=="strongly disagree"] <- 1

dta3$BLACKPOLICIES_1S <- NA
dta3$BLACKPOLICIES_1S[dta3$blackpolicies_1=="Strongly Agree"] <- 4
dta3$BLACKPOLICIES_1S[dta3$blackpolicies_1=="Somewhat Agree"] <- 3
dta3$BLACKPOLICIES_1S[dta3$blackpolicies_1=="Somewhat Disagree"] <- 2
dta3$BLACKPOLICIES_1S[dta3$blackpolicies_1=="Strongly Disagree"] <- 1

dta3$BLACKPOLICIES_7S <- NA
dta3$BLACKPOLICIES_7S[dta3$blackpolicies_7=="strongly agree"] <- 4
dta3$BLACKPOLICIES_7S[dta3$blackpolicies_7=="somewhat agree"] <- 3
dta3$BLACKPOLICIES_7S[dta3$blackpolicies_7=="somewhat disagree"] <- 2
dta3$BLACKPOLICIES_7S[dta3$blackpolicies_7=="strongly disagree"] <- 1

dta3$BLACKPOLICIES_8S <- NA
dta3$BLACKPOLICIES_8S[dta3$blackpolicies_8=="strongly agree"] <- 4
dta3$BLACKPOLICIES_8S[dta3$blackpolicies_8=="somewhat agree"] <- 3
dta3$BLACKPOLICIES_8S[dta3$blackpolicies_8=="somewhat disagree"] <- 2
dta3$BLACKPOLICIES_8S[dta3$blackpolicies_8=="strongly disagree"] <- 1

dta3$BLACKFAVORS_1S <- NA
dta3$BLACKFAVORS_1S[dta3$blackfavors_1=="Strongly Agree"] <- 4
dta3$BLACKFAVORS_1S[dta3$blackfavors_1=="Somewhat Agree"] <- 3
dta3$BLACKFAVORS_1S[dta3$blackfavors_1=="Somewhat Disagree"] <- 2
dta3$BLACKFAVORS_1S[dta3$blackfavors_1=="Strongly Disagree"] <- 1

dta3$BLACKFAVORS_7S <- NA
dta3$BLACKFAVORS_7S[dta3$blackfavors_7=="strongly agree"] <- 4
dta3$BLACKFAVORS_7S[dta3$blackfavors_7=="somewhat agree"] <- 3
dta3$BLACKFAVORS_7S[dta3$blackfavors_7=="somewhat disagree"] <- 2
dta3$BLACKFAVORS_7S[dta3$blackfavors_7=="strongly disagree"] <- 1

dta3$BLACKFAVORS_8S <- NA
dta3$BLACKFAVORS_8S[dta3$blackfavors_8=="strongly agree"] <- 4
dta3$BLACKFAVORS_8S[dta3$blackfavors_8=="somewhat agree"] <- 3
dta3$BLACKFAVORS_8S[dta3$blackfavors_8=="somewhat disagree"] <- 2
dta3$BLACKFAVORS_8S[dta3$blackfavors_8=="strongly disagree"] <- 1

dta3$CTBLKHS <- NA
dta3$CTBLKHS[dta3$HS3_8=="Never or almost never"] <- 0
dta3$CTBLKHS[dta3$HS3_8=="Less than once a month"] <- 1/60
dta3$CTBLKHS[dta3$HS3_8=="1-3 times each month"] <- 2/30
dta3$CTBLKHS[dta3$HS3_8=="At least once a week"] <- 4/30
dta3$CTBLKHS[dta3$HS3_8=="Every day"] <- 1
dta3$CTBLKHS[dta3$HS3_8=="Refused"] <- NA

dta3$CTBLKHSL <- NA
dta3$CTBLKHSL[dta3$HS3_8=="Never or almost never"] <- 0
dta3$CTBLKHSL[dta3$HS3_8=="Less than once a month"] <- 1
dta3$CTBLKHSL[dta3$HS3_8=="1-3 times each month"] <- 2
dta3$CTBLKHSL[dta3$HS3_8=="At least once a week"] <- 3
dta3$CTBLKHSL[dta3$HS3_8=="Every day"] <- 4
dta3$CTBLKHSL[dta3$HS3_8=="Refused"] <- NA

dta3$SEEBLK7 <- NA
dta3$SEEBLK7[dta3$RE28_7=="Never or almost never"] <- 0
dta3$SEEBLK7[dta3$RE28_7=="Less than once a month"] <- 1/60
dta3$SEEBLK7[dta3$RE28_7=="1-3 times each month"] <- 2/30
dta3$SEEBLK7[dta3$RE28_7=="At least once a week"] <- 4/30
dta3$SEEBLK7[dta3$RE28_7=="Every day"] <- 1
dta3$SEEBLK7[dta3$RE28_7=="Refused"] <- NA

dta3$TLKBLK7 <- NA
#dta3$SEEBLK7[dta3$RE29_7==""]
dta3$TLKBLK7[dta3$RE29_7=="Never or almost never"] <- 0
dta3$TLKBLK7[dta3$RE29_7=="Less than once a month"] <- 1/60
dta3$TLKBLK7[dta3$RE29_7=="1-3 times each month"] <- 2/30
dta3$TLKBLK7[dta3$RE29_7=="At least once a week"] <- 4/30
dta3$TLKBLK7[dta3$RE29_7=="Every day"] <- 1
dta3$TLKBLK7[dta3$RE29_7=="Refused"] <- NA

dta3$TLKBLK7L <- NA
#dta3$SEEBLK7[dta3$RE29_7==""]
dta3$TLKBLK7L[dta3$RE29_7=="Never or almost never"] <- 0
dta3$TLKBLK7L[dta3$RE29_7=="Less than once a month"] <- 1
dta3$TLKBLK7L[dta3$RE29_7=="1-3 times each month"] <- 2
dta3$TLKBLK7L[dta3$RE29_7=="At least once a week"] <- 3
dta3$TLKBLK7L[dta3$RE29_7=="Every day"] <- 4
dta3$TLKBLK7L[dta3$RE29_7=="Refused"] <- NA

dta3$HELPBLKS_7 <- NA
dta3$HELPBLKS_7[dta3$RE23_2_7=="Government should help blacks1"] <- 7
dta3$HELPBLKS_7[dta3$RE23_2_7=="2"] <- 6
dta3$HELPBLKS_7[dta3$RE23_2_7=="3"] <- 5
dta3$HELPBLKS_7[dta3$RE23_2_7=="4"] <- 4
dta3$HELPBLKS_7[dta3$RE23_2_7=="5"] <- 3
dta3$HELPBLKS_7[dta3$RE23_2_7=="6"] <- 2
dta3$HELPBLKS_7[dta3$RE23_2_7=="Blacks should help themselves7"] <- 1

dta3$HELPBLKS_8 <- NA
dta3$HELPBLKS_8[dta3$RE23_2_8=="Government should help blacks1"] <- 7
dta3$HELPBLKS_8[dta3$RE23_2_8=="2"] <- 6
dta3$HELPBLKS_8[dta3$RE23_2_8=="3"] <- 5
dta3$HELPBLKS_8[dta3$RE23_2_8=="4"] <- 4
dta3$HELPBLKS_8[dta3$RE23_2_8=="5"] <- 3
dta3$HELPBLKS_8[dta3$RE23_2_8=="6"] <- 2
dta3$HELPBLKS_8[dta3$RE23_2_8=="Blacks should help themselves7"] <- 1

dta3$SOUTH8 <- 1*(dta3$PPREG4_8 %in% ("South"))
dta3$MARRIED7 <- 1*(dta3$PPMARIT_7=="Married")
dta3$METRO7 <- 1*(dta3$PPMSACAT_7=="Metro")

dta3$RETIRED7 <- 1*(dta3$PPWORK_7=="Not working - retired")
dta3$DISABLED7 <- 1*(dta3$PPWORK_7=="Not working - disabled")

dta3$STCOFIPSRE <- dta3$STFIPSRE*1000+dta3$COFIPSRE

dta3$STAY <- NA
dta3$STAY[dta3$STCOFIPSRE==dta3$STCOFIPS & ! dta3$STCOFIPSRE %in% c(NA)] <- 1
dta3$STAY[(! dta3$STCOFIPSRE==dta3$STCOFIPS) & ! dta3$STCOFIPSRE %in% c(NA)] <- 0

dta3$OBAMA5 <- 1*(dta3$VP9_5=="Barack Obama and Joe Biden, the Democrats")
dta3$OBAMA3 <- 1*(dta3$VP5_3=="Barack Obama, the Democrat")
dta3$OBAMA4 <- 1*(dta3$VP5_4==2)

dta3$COUNEMP10 <- dta3$SE_T033_006/(dta3$SE_T033_005+dta3$SE_T033_006)

dta3$BLKSCALE87 <- dta3$BLKSCALE8-dta3$BLKSCALE7
dta3$BLKSCALE87B <- (dta3$BLKSCALE87 > 0)*1

dta3$BLKSCALE83 <- dta3$BLKSCALE8-dta3$BLKSCALE3
dta3$BLKSCALE83B <- (dta3$BLKSCALE83 > 0)*1

dta3$BLKSCALE84 <- dta3$BLKSCALE8-dta3$BLKSCALE4
dta3$BLKSCALE84B <- (dta3$BLKSCALE84 > 0)*1


####### PREJUDICE MEASURE WITH ESTEEM
dta3$DETH1_A_3B <- dta3$DETH1_A_3
dta3$DETH1_A_3B[dta3$DETH1_A_3 < 0] <- NA

dta3$DETH1_B_3B <- dta3$DETH1_B_3
dta3$DETH1_B_3B[dta3$DETH1_B_3 < 0] <- NA

dta3$DETH1_C_3B <- dta3$DETH1_C_3
dta3$DETH1_C_3B[dta3$DETH1_C_3 < 0] <- NA

dta3$WHTHRDWRK3 <- as.numeric(dta3$DETH1_A_3B)
dta3$WHTINTEL3 <- as.numeric(dta3$DETH1_B_3B)
dta3$WHTTRST3 <- as.numeric(dta3$DETH1_C_3B)
dta3$WHTSCALE3 <-  dta3$WHTTRST3+dta3$WHTINTEL3+dta3$WHTHRDWRK3

dta3$WHTHRDWRK4 <- as.numeric(dta3$DETH1_A_4)
dta3$WHTINTEL4 <- as.numeric(dta3$DETH1_B_4)
dta3$WHTTRST4 <- as.numeric(dta3$DETH1_C_4)
dta3$WHTSCALE4 <-  dta3$WHTTRST4+dta3$WHTINTEL4+dta3$WHTHRDWRK4

dta3$WHTHRDWRK5 <- as.numeric(dta3$DETH1_A_5)
dta3$WHTINTEL5 <- as.numeric(dta3$DETH1_B_5)
dta3$WHTTRST5 <- as.numeric(dta3$DETH1_C_5)
dta3$WHTSCALE5 <-  dta3$WHTTRST5+dta3$WHTINTEL5+dta3$WHTHRDWRK5

dta3$WHTHRDWRK6 <- as.numeric(dta3$DETH1_A_6)
dta3$WHTINTEL6 <- as.numeric(dta3$DETH1_B_6)
dta3$WHTTRST6 <- as.numeric(dta3$DETH1_C_6)
dta3$WHTSCALE6 <-  dta3$WHTTRST6+dta3$WHTINTEL6+dta3$WHTHRDWRK6

dta3$WHTHRDWRK7 <- as.numeric(dta3$DETH1_A_7)
dta3$WHTINTEL7 <- as.numeric(dta3$DETH1_B_7)
dta3$WHTTRST7 <- as.numeric(dta3$DETH1_C_7)
dta3$WHTSCALE7 <-  dta3$WHTTRST7+dta3$WHTINTEL7+dta3$WHTHRDWRK7

dta3$WHTHRDWRK8 <- as.numeric(dta3$DETH1_A_8)
dta3$WHTINTEL8 <- as.numeric(dta3$DETH1_B_8)
dta3$WHTTRST8 <- as.numeric(dta3$DETH1_C_8)
dta3$WHTSCALE8 <-  dta3$WHTTRST8+dta3$WHTINTEL8+dta3$WHTHRDWRK8


dta3$BLKHRDWRK5 <- NA
dta3$BLKHRDWRK5[! dta3$PPETHM_5=="Black, Non-Hispanic" & ! dta3$PPETHM_5 %in% c(NA)] <-as.numeric(dta3$DETH2_A_5[! dta3$PPETHM_5=="Black, Non-Hispanic" & ! dta3$PPETHM_5 %in% c(NA)])
dta3$BLKHRDWRK5[dta3$DETH2_A_5 < 0] <- NA

dta3$BLKINTEL5 <- NA
dta3$BLKINTEL5[! dta3$PPETHM_5=="Black, Non-Hispanic" & ! dta3$PPETHM_5 %in% c(NA)] <-as.numeric(dta3$DETH2_B_5[! dta3$PPETHM_5=="Black, Non-Hispanic" & ! dta3$PPETHM_5 %in% c(NA)])
dta3$BLKINTEL5[dta3$DETH2_B_5 < 0] <- NA

dta3$BLKTRST5 <- NA
dta3$BLKTRST5[! dta3$PPETHM_5=="Black, Non-Hispanic" & ! dta3$PPETHM_5 %in% c(NA)] <-as.numeric(dta3$DETH2_C_5[! dta3$PPETHM_5=="Black, Non-Hispanic" & ! dta3$PPETHM_5 %in% c(NA)])
dta3$BLKTRST5[dta3$DETH_C_5 < 0] <- NA

dta3$BLKSCALE5 <- dta3$BLKHRDWRK5+dta3$BLKINTEL5+dta3$BLKTRST5


dta3$BLKHRDWRK6 <- NA
dta3$BLKHRDWRK6[! dta3$PPETHM_6=="Black, Non-Hispanic" & ! dta3$PPETHM_6 %in% c(NA)] <-as.numeric(dta3$DETH2_A_6[! dta3$PPETHM_6=="Black, Non-Hispanic" & ! dta3$PPETHM_6 %in% c(NA)])
dta3$BLKHRDWRK6[dta3$DETH2_A_6 < 0] <- NA

dta3$BLKINTEL6 <- NA
dta3$BLKINTEL6[! dta3$PPETHM_6=="Black, Non-Hispanic" & ! dta3$PPETHM_6 %in% c(NA)] <-as.numeric(dta3$DETH2_B_6[! dta3$PPETHM_6=="Black, Non-Hispanic" & ! dta3$PPETHM_6 %in% c(NA)])
dta3$BLKINTEL6[dta3$DETH2_B_6 < 0] <- NA

dta3$BLKTRST6 <- NA
dta3$BLKTRST6[! dta3$PPETHM_6=="Black, Non-Hispanic" & ! dta3$PPETHM_6 %in% c(NA)] <-as.numeric(dta3$DETH2_C_6[! dta3$PPETHM_6=="Black, Non-Hispanic" & ! dta3$PPETHM_6 %in% c(NA)])
dta3$BLKTRST6[dta3$DETH_C_6 < 0] <- NA

dta3$BLKSCALE6 <- dta3$BLKHRDWRK6+dta3$BLKINTEL6+dta3$BLKTRST6

dta3$PREJUDICE3B <- dta3$WHTSCALE3-dta3$BLKSCALE3
dta3$PREJUDICE4B <- dta3$WHTSCALE4-dta3$BLKSCALE4
dta3$PREJUDICE5B <- dta3$WHTSCALE5-dta3$BLKSCALE5
dta3$PREJUDICE6B <- dta3$WHTSCALE6-dta3$BLKSCALE6
dta3$PREJUDICE7B <- dta3$WHTSCALE7-dta3$BLKSCALE7
dta3$PREJUDICE8B <- dta3$WHTSCALE8-dta3$BLKSCALE8
#####
